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Abstract 

Dark matter modeled as a classical scalar field that interacts only with gravity and with itself by a potential that is 
close to quartic at large field values and approaches a quadratic form when the field is small would be gravitationally 
produced by infiation and at the present epoch could act like an ideal fluid with pressure that is a function only 
of the mass density. This could have observationally interesting effects on the core radii and solid body rotation of 
dark matter halos and on the low mass end of the primeval mass fluctuation power spectrum. 
Subject headings: cosmology: dark matter — galaxies: structure 



1. Introduction 

This paper continues the discussion by Peebles & 
Vilenkin (1999, hereafter PV) of a model for dark 
matter as a gravitationally produced self-interacting 
scalar field. It is shown that this dark matter acts like 
an ideal fluid with pressure that is a function of den- 
sity alone, with possibly interesting consequences for 
the structure of dark halos of galaxies at low redshift 
and the evolution of the dark matter distribution at 
redshift z 10''^ out of which the galaxies are thought 
to have formed. These effects are discussed in §2, and 
the fluid model derived in §3. The distinction between 
the behavior of this fluid and a gas of interacting dark 
matter particles (Spergel & Steinhardt 1999) is dis- 
cussed in §4, along with the possible significance of 
the fluid picture. 

The model is motivated by the following consider- 
ations. Fields that interact only with gravity would 
be gravitationally produced by inflation (Ford 1987), 
and a gravitationally produced scalar fleld would be 
left in a squeezed state (Grishchuk & Sidorov 1990) 
that acts like a classical field. One could imagine that 
during inflation the potential of a scalar field that will 
end up as dark matter is close toV = and that 

inflation is driven by a fleld with another quartic 
potential, U = A0(/)'*/4. If the dimensionless coefla- 
cients satisfy (PV eqs. [85] and [86]; for discussions of 
other aspects of this constraint see Kofman & Linde 
1987; Felder, Kofman & Linde 1999; and Starobinsky 
& Yokoyama 1994) 
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the mass density py in dark matter at the end of infla- 
tion may be an interesting fraction of the mass density 
that becomes ordinary matter and radiation, and the 
space distribution of py is smooth enough that it does 



not violate the isotropy of the thermal background ra- 
diation. When V{y) is close to quartic the dark mass 
density after inflation varies as py oc a{t)~'^, where 
a{t) is the expansion parameter, meaning py is a fixed 
fraction of the mass density in radiation. Since the 
field y is supposed to have low pressure now the curva- 
ture of the potential has to become close to quadratic 
at small \y\. I adopt the form 



V = m^yy2 + K\y\yq, q>2, 



(2) 



for typical field values at redshifts z < 10^. PV take 
5 = 4, making K dimensionless. The extra parameter 
q broadens the possibilities for fitting to the astron- 
omy. 

The analysis in |3 shows that when the anharmonic 
term in equation (0) is subdominant the y-field mass 
behaves like an ideal nonrelativistic fluid with density 
Py = m?{y'^) and pressure py = Fpy''/'^ , where 



F 



q 



2 29/2-iXr((g+l)/2) 



j,l/2^q T{{q + 2)/2)' 



(3) 



The velocity of sound, Cg = {dpy/dpyY/^, multiplied 
by the gravitational time (Gpy)^^^"^ , is a Jeans length 

9/4-1 

(X Py' 

2. The Jeans Length 

In a pressure-supported dark matter halo the mass 
density as a function of radius r near the center is 



Pyir) = Pc{l-r^/rl), = 



3ci 



(4) 



27rGp, 

The rotation curve at r < is the solid body form 

v{r)^Vcr/rc, Uc = (5) 
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Since the halo can be supported by Reynolds stress 
as well as pressure the condition on the size of the 
region of solid body rotation in a halo of fluid dark 
matter with central density pc is 



rc > Rp\ 
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(6) 



for some constant, R. 

A second constraint follows from the evolution 
of the fluctuation spectrum of the j/-mass distribu- 
tion through the end of radiation-dominated expan- 
sion at redshift Zcq = 2.4 x XQ'^^yh? and world time 
teq ~ 3 X lQ^^Q.~'^h~'^ s, for density parameter fi^ 
in dark matter and Hubble parameter h in units of 
100 km s~^ Mpc""'^. The comoving wavenumber kc 
(scaled to physical length at the present epoch) at 
the first zero of the transfer function produced by the 
y-field pressure satisfies, in order of magnitude (Pee- 
bles 1980, eq. [92.27]), fccCstoq^oq ~ tt. At fc ^ fcc the 
transfer function is suppressed relative to pressureless 
dark matter. A lower bound on kc for an acceptable 
model of structure formation yields an upper bound 
on at z ~ Zcq- 

At the bounds on q from kc and equation (^, 



„„„„riy/i^ 50 km s 
3000- ^ 



kcVc 1 Mpc 



(7) 



3 X IQ'^^lK 



g / hvc 50 km s ^ 
\vc 1 kpc 



21 (9-2)/4 



From the simulations in White & Kroft (2000) it ap- 
pears that if kc ^ lO/i Mpc~^ the transfer function 
might be acceptable for structure formation. At this 
bound, with ily = 0.25 and h = 0.7, and halo pa- 
rameters Vc = 50 km s~^ at Vc = \h^^ kpc, equa- 
tion (|^) says q — 3.7. Increasing the Jeans length to 
hrc = 3 kpc gives q = 3.4. For Oy = 1 with the other 
parameters from the first example, q = 3.4. If g = 3.6 
the physical Jeans length varies quite slowly with the 
mass density, Tc oc Py^'^- liq = A the Jeans length is a 
universal constant: hvc <, 0.5 kpc at kc ^ lO/i Mpc"""^, 
Vty = 0.25, and h = 0.7. 

3. Fluid Dark Matter 

The method of analysis for a pure quartic poten- 
tial, and numerical tests of the approximations, are 
presented in Peebles (1999). The present case as- 
sumes the velocity of sound and the dark matter 
streaming velocity are nonrelativistic. The former im- 
plies the anharmonic part of V{y) is subdominant, so 



the field is oscillating at frequency m. li q ^ A the 
equation of state is Py ~ iXp^ / {^rn^) , where K ~ X 
(eq, §). and 
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(8) 



It seems safe to assume this is much more rapid than 
the rate at which the mass distribution and flux are 
changing, so we can find the mass density and ve- 
locity as functions of position and time by averaging 
functions of the field over a time interval r such that 



Py/Py < 



^ m. 



(9) 



Thus the mass density is well approximated as 

py^{f}^m^{y'). (10) 

In the time and space intervals t and vt a good ap- 
proximation to a solution to the wave equation is 



y = 77(7(t-vr)). 



(11) 



where = —dV{ri)/dri. This is a Lorentz transforma- 
tion of ri{t), with Lorentz factor 7, so it describes 
proper mass density (77^) moving with velocity v. 
Consistent with this, the field gradient is 

Vy = -yv, (12) 

for small v, and the mass flux density is 

{iy) = -{y\/y)^Py^. (13) 

Consider first the pressure balance in a static dark 
matter halo. The line element is ds^ = (1 + 2(j))dt'^ — 
(1 — 2(j>)dr^, where the Newtonian potential (p is a. 
function of position but independent of time, and 
terms of order (f)'^ dropped. The wave equation is 

(1 - A^)y - V^y 

+ (l-20)(m2y + Xy|yr-2)=O. (14) 

The result of multiplying this equation by y and av- 
eraging over time, and using {yy) = —{if'), from in- 
tegration by parts, is the virial relation 

{\-A4>){f) + {yV^y) 

= (l-2</.)(m2(y2) + if(|y|?)). (15) 

The result of multiplying equation ( [T^ by the gradi- 
ent of y with respect to the Cartesian position compo- 
nent Tq, averaging over time, and using equation (Eq), 
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IS 



.„(m2y2 + X|y|''-2j/V2j/) 
+a;3((2/,.2/,^)-9„5^(y'>/4) 
+ (1 - 2(b)d^py = 0, 



(16) 



where 

Py = iq-2)K{\yn/{2q). (17) 

The coefBcient of (f>^a may be written as w?{y^) = py, 
because the anharmonic term is subdominant in a 
nonrelativistic fluid and equations ( p^ ) and (12) say 
(yV^y) ^ {yyv^) ^ —pyV^. The second term in the 
second line is small because it depends on the astro- 
nomical length scale of variation of the mean density. 
In last term the Newtonian potential is an unimpor- 
tant correction to the pressure gradient force. Thus 
the balance condition is 



Py(j)^a + daPy + d/sipyVaVp) = 0. 



(18) 



The last term, which uses equation (^^, is the ana- 
log of the Reynolds stress in a turbulent fluid or the 
velocity dispersion in a gas. Equation ( p^ ) for py is 
reduced to equation (||) by using y oc sin mi. 

When the Reynolds stress term may be neglected 
equation ( p!8| ) yields the familiar relation between the 
halo central density, the core radius Tc, and the ve- 
locity of sound in equation (||) . 

Consider next dynamics under the simplifying as- 
sumption that the time-dependence of the Newtonian 
potential (j) may be ignored. The virial relation in 
equation (^5|) applies here as an average over the time 
and space intervals t and vt (eq. |^). The result of 
using the wave equation to eliminate y from the time 
derivative of the mass flux density f = — f^f, aver- 
aging, and using the virial relation, is 

d{fc,)/dt + <j>A^^y^ + K\yV - 2yV2j/) 
+ {l + A4>)dp{{y,o.y,p)-dMy^)/A) 
+ (1 + 2^)dc.py = 0. (19) 

This generalizes equation (|l^ ) to a dynamical situa- 
tion. As before, the nonrelativistic limit is 

d{pyVa)/dt + dpipyVaVp) = -Py4>,a - OaPy (20) 

The mass density is 

- (1 - 40)yV2 + VyVS 
+ {1^2^){m^yy2 + K\y\yq), (21) 



with Py = —V • f = —V • pyV. This with equation (20) 
is the Euler equation — the Navier-Stokes equation 
with zero viscosity — for a fluid with pressure Py in 
the gravitational field g = — V0. 

Similar methods show the mean mass density in an 
expanding universe with ^ = satisfies 



(22) 



Here again the intrinsic pressure (eq. [|T^) is added 
to the contribution {pyv'^)/3 by turbulence in a fluid 
or velocity dispersion in a gas. 

The ideal fluid picture cannot be complete: the 
velocity of sound is an increasing function of py, so 
pressure waves tend to grow into the analog of shock 
waves. The numerical examples in Peebles (1999) 
for a pure quartic potential indicate that, unlike real 
shock waves, the shock-like features tend to be re- 
versible. This fluid picture with reversible "shocks" 
must eventually break down by the cascade of field 
energy to wavelengths short enough that the better 
analog is a gas of particles. To be investigated is the 
time for this to happen. This limitation of the fluid 
picture does not seem to be a problem for the halo 
model (eq. iQ), because it depends on the apparently 
reasonable assumption that the y-field is a stationary 
time process. 

4. Discussion 

Dark matter pictured as an interacting gas has 
been considered by Carlson, Machacek, & Hall (1992), 
Machacek (1994), and de Laix, Scherrer, & Schaefer 
(1995). Spergel & Steinhardt (1999) show interac- 
tions can have interesting effects on the dark matter 
distribution; the remarkable burst of papers on the 
concept includes Ostriker (1999); Hannestad (1999); 
Miralda-Escude (2000); Moore et al. (2000); Hogan 
& Dalcanton (2000); Yoshida et al. (2000); Burkert 
(2000); and Firmani et al. (2000). This gas picture 
assumes small quantum occupation numbers. The 
fiuid picture discussed here and in PV assumes the 
opposite limit of dark matter squeezed — perhaps by 
inflation — into a state with large occupation num- 
bers that behaves like a classical field. The quadratic 
part of V{y) (eq. [^) describes particles of mass m 
and present mean number density lOOA^^/'* cm^^ 
(eq. [|j), but one should think of these particles as 
moving coherently with large de Broglie wavelength, 
as a classical field. The examples in §3 show this field 
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can behave like an ideal fluid. 

The interacting gas picture may be approximated 
as a fluid, too, but with two important differences. 
First, the Jeans length in a patch of self-interacting 
gas depends on its history, while the potential V{y) 
fixes an intrinsic Jeans length (eq. |^). Second, self- 
interacting gas acts as a viscous fluid, with ram pres- 
sure, while the analysis here and in Peebles (1999) 
shows situations where the fluid dark matter has no 
viscosity (eq. [Q), and soliton-likc behavior (Peebles 
1999). 

Pressure in the fluid dark matter could suppress 
the small-scale end of the spectrum of primeval fluctu- 
ations in the dark mass distribution. Kamionkowski 
& Liddle (1999) suggest a suppression could correct 
indications of excess numbers of dwarf galaxies in 
the adiabatic cold dark matter model (Klypin et al. 
1999 and references therein), and note the suppres- 
sion could have originated during inflation. Warm 
dark matter (Hogan & Dalcanton 2000) and fluid dark 
matter produce a similar effect from a simpler infla- 
tion, at the price of more complicated dark matter. 

The model for the potential (eq. 1^) is most ele- 
gant if the anharmonic part is quartic, q = 4, so V{y) 
is an analytic function of y that can directly relate to 
a simple model for eternal inflation (eq. ||l|). In this 
case natural choices for the other parameters require 
Jeans length hrc ^ 0.5 kpc, which may be somewhat 
small for dark matter halos. Halos can be supported 
at larger core radii in the ordinary way, by Reynolds 
stress (eq. Q ) , but this leaves open the issue of solid 
body rotation curves. Since V{y) is an effective po- 
tential for the classical field limit it might not be un- 
reasonable to consider non-integer q. IfO<g — 2<c2 
then V{y) can produce dark matter halos with inter- 
esting solid body rotation without seriously affecting 
the primeval mass fluctuation spectrum. If g ~ 3.5 
the model can have interesting effects on both. 

The fluid model requires the core radii Vc of dark 
halos satisfy equation (^ as a lower bound; veloc- 
ity dispersion can only increase Vc at given central 
density. Dalcanton & Bernstein (1999) present an el- 
egant example of solid body rotation in a low surface 
brightness galaxy that seems to be a good approxi- 
mation to a dark matter halo. The critical question, 
whether such behavior is universal, has a long history 
and is under discussion (Salucci & Persic 1999; van 
den Bosch et al. 1999; Swaters, Madore, & Trehwella 
2000; Firmani et al. 2000; and references therein). 



The survival of the dark matter in dwarf satellite 
galaxies is an important issue for the interacting gas 
picture (Moore et al. 2000) and the fluid picture. In 
the simple cases in Peebles (1999) lumps in the fluid 
dark matter mimic solitons. Fluid dark matter in the 
halos of dwarf satellite galaxies might similarly pass 
through a larger dark matter halo with little distur- 
bance, but the idea has not been considered in any 
detail. 

One may also judge models by elegance. Thus 
one might seek a simple reason why the present ratio 
Pin/ Pr of mass densities in dark matter and radiation 
is not far from unity, compared to the great difference 
at the end of inflation if the dark matter were quite 
cold and noninteracting. This also has a long history, 
as in Lee & Weinberg (1977); a recent proposal by 
Zlatev & Steinhardt (1999) is another model of inter- 
acting dark matter, with potential V{y) such that the 
present value of Pm/ Pr is insensitive to a broad range 
of initial conditions. The potential in equation yj^ 
does not exhibit this behavior but does allow another 
consideration. It is easy to imagine the field that 
acted as the inflaton during inflation had energy den- 
sity only a few orders of magnitude larger than some 
other field (Kofman & Linde 1987), and that the lat- 
ter field happens to interact only with itself and grav- 
ity, so it ends up as dark matter (PV). The quartic 
form of V{y) at large field values would keep Pm/ Pr 
close to unity. The potential must change to near 
quadratic at the rms value of y reached some twenty 
orders of magnitude of expansion after infiation ends, 
in time to make the dark mass density comparable 
to the baryons. This still is a considerable "cosmic 
coincidence." 

The bound on A (eq. Q) at high redshift depends 
on the highly uncertain efficiency of production of or- 
dinary matter and radiation as infiation ends. Within 
this uncertainty A may be adjusted so structure for- 
mation is seeded by the curvature fluctuations from 
the inflaton, as in the adiabatic cold dark matter 
model. The replacement of cold with fluid dark mat- 
ter in this model would affect structure formation only 
on relatively small scales, where it might be benefi- 
cial. Alternatively, one may postulate a larger value 
of A that would produce significant isocurvature fluc- 
tuations, perhaps in a more elaborate version of the 
model in Hu & Peebles (1999). This could be ben- 
eflcial if the model assembled galaxies earlier than 
predicted in the cold dark matter model, and earlier 
formation better fit the observations. 
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